Magnetized Kelvin-Helmholtz instability in the presence 
of a radiation field 
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Abstract The purpose of this study is to analyze the 
dynamical role of a radiation field on the growth rate 
of the unstable Kelvin - Hclmholtz (KH) perturbations. 
As a first step toward this purpose, the analyze is done 
in a general way, irrespective of applying the model 
to a specific astronomical system. The transition zone 
between the two layers of the fluid is ignored. Then, 
we perform a linear analysis and by imposing suitable 
boundary conditions and considering a radiation field, 
we obtain appropriate dispersion relation. Unstable 
modes are studied by solving the dispersion equation 
numerically, and then growth rates of them are ob- 
tained. By analyzing our dispersion relation, we show 
that for a wide range of the input parameters, the radi- 
ation field has a destabilizing effect on KH instability. 
In eruptions of the galaxies or supermassive stars, the 
radiation field is dynamically important and because of 
the enhanced KH growth rates in the presence of the 
radiation; these eruptions can inject more momentum 
and energy into their environment and excite more tur- 
bulent motions. 

Keywords Kelvin-Helmholtz instability; radiation; 
compressible media; magnetohydrodynamic. 



1 Introduction 

Kelvin-Helmholtz (KH) instability occurs when dif- 
ferent layers of a fluid are in relative motion. This 
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instability occurs in many astrophysical systems and 
has bee n widely stud ied in relation to the solar 



wind (lAmerstorfer et al.l 120071 : iBettarini et al.l 12 006 



Hasegawa et al.ll2004 1. puls ar winds (Bucciantini and Del Zannafr 



2006t ) I and thermal flares (jVenter and Meintiesl 120061 ) 
The KH instability can also drive mixing and tur- 
bulence and is thus r elevant in protoplanetary discs 
( Johansen et al. 20061 Gomez and Ostriker 20051). ac- 
cretion discs and magnetosphere s (ILi and Naravan 



2004 ) and other jets and outflows (jBatv and Keppens 
2006 L 



Chandrasekharl (|196lh studied KH instability in the 
linear regime. Using perturbation analysis, he de- 
rived the conditions and the timescales required for 
the growth of this instability. Also examined were 
the effects of other parameters such as radiative cool- 
ing and magnetic fields, and their ability to stabi- 
lize the development of these non-linear structures. 
More recently, KH instability has been investigated by 
means of numerical simulations, in both the relativistic 



( Perucho et al. 2004 , 20061 Bucciantini and Del Zanna 


2006f) and non-relativistic limits ( 


Rossi et al. 1997; 


Downes and Ravll998HMicono et al. 


1998; Hanasz and Soli 


19981 Micono et al. 200(J Michikoshi 


and Inutsukal2006l).| 



One very important factor in all KH simulations is the 
role played by the magnetic field. Magnetohydrody- 
namics is required for all numerical simulations in or- 
der to accurately repre s ent the presence o f this field 



( Viallet and Batvl l2007t IJeong et al 



2000; 



Min 1997 



Shadmehri and Downesl2008t Shadmehri and Rammos 
2010h . 

But radiation may have dynamical effects in some 
of the astronomical objects. For example, formation of 
the stars takes place in the interstellar medium (ISM) 
where a wide variety of the structures and the physi- 
cal conditions are observed. There is a growing interest 
toward understanding properties of the ISM, in partic- 
ular regarding to the dynamical role of the radiation in 
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the gravit ational instabi l ity and possibly massive star 
formation ( Vraniesl Il99d: IVranies and Cadezlll990l) . 

Dynamical role of the radiation field in the context 
of the accretion flows has also been studied by many 
authors. At the larger scales, we have radiation domi- 
nated discs. In particular, the very inner parts of the ac- 
cretion discs around black holes are radiation pressure 
dominated, though there are many uncertainty about 
their p hysical state. In the standard model of accretio n 



discs (jBath et al.l 11974 IShakura and Sunvaevl 119731 ). 



the viscous stress is assumed to scale with the radiation 
pressure, but the disc itsel f is subject to the thermal and 
the viscous instabilities (IShakura and Sunvaevl 11976 ; 
Lightman and Eardlevl 119741) . In addition to these 
instabilities, dynamical instabilities also exist. First 
and foremost, magnetor otational instability (MRI) 



the standard linear approach of KH instability to in- 
clude heating and cooling. However, dynamical role of 
the radiation field has not been studied analytically or 
numerically as far as wc know. Considering the impor- 
tance of the radiation field in some of the astronomical 
systems subject to the KH instability as we described 
above, our goal is to investigate KH instability includ- 
ing the dynamical role of the radiation. In the next 
section, basic assumptions and equations are presented. 
In section 3, wc will study KH instability in the linear 
regime using a perturbation analysis. Different configu- 
rations are considered and we will generalize our study 
to include not only the radiation field but the magnetic 
field. 



(jHawlev and Balbua 119911) is effective in weakly mag- 2 Basic Equations 



netized accretion discs. Gammie (1998) has suggested 
that th e overstable photon bubble modes discussed by 
AronsJ(1992) in the X-ray pulsars also exist in radi- 
ation pressure dominated accretion flows in general. 
These dynamical instabilities may all play a role at 
some level in the dynamics and thermodynamics of the 
radiation pressure domina ted portion of accretion disks 
( Blaes and Socratesll200ll ). 

According to the observations of the galaxies only 
a few percent of the available gas reservoir galax- 
ies is converted into st a rs per local fre e-fall time 
( Krumholz and Tanl [20071 iKennicutt] 1 19981 ). In other 



words, the efficiency of star formation in the galaxies is 
low. Injection of energy and momentum into the inter- 
stellar medium (ISM) by stellar processes (feedback), 
may be responsible for the inefficiency of star formation 
in galaxies and for their self-regulation. Various models 
for feedback, effects of HII gas pressure, shocked stellar 
winds, protostellar jets have been proposed over re- 
cent years. Dust grains are hydrodynamically coupled 
to the gas, and then radiation pressure on dust can 
help stabilize the gas against its own self-gravity and 
may therefore be an i mportant feedback mechanism 
( Thompson et al J 12005). When a jet or outflow inter- 
acts with its ambient medium, not only mass but energy 
are exchanged via the induced turbulence. One of the 
mechanisms that may initiate such a turbulence is KH 
instability. It has been proposed that winds may launch 
via a radiation-driven m echanisms (jSharma and Nath 
201 it iMurrav et aHl2010| ). Considering the dynamical 
effect of radiation in such system, it would be interest- 
ing to study the resulting KH instability at its interface 
with the ambient medium. 

Early analytical studies of KH instability generally 
neglect an imbalance of the cooling and heating mech- 
anisms. But subsequent analytical studies generalized 



Our main equations are the standard MHD equa- 
tions including a radiation field assuming that the gas 
and the radiation are i n per fe ct radiative equilibrium 
( Mihalas and Mihalasl 1984 ; Shadmehri and Downed 
20071). Thus" 



^ + V-(pu) = 0, 



^+U'Vu = --V( Pg +p r ) + —(J x B), 



/9R 

- = Vx(uxB), 



V-B = 0, 



dE 
~dt 



u-VE+(E + P)V • u = -V • F, 



F = -— Ve r , 

6Kp 



Pg = P C s 



(1) 

(2) 
(3) 
(4) 

(5) 
(6) 
(7) 



where the subscripts r and g refer to the radiation and 
the gas components, respectively. Here, p and u are the 
density and the velocity of the gas. Also, the gas and 
the radiation and the total pressures are denoted by p s 
and p r and P = p g + p r , respectively. Moreover, we 
introduce the total internal energy E = e g + e r , where 
e g = Pg/i'J— 1) and e r = 3p r and 7 is the adiabatic index 
of the gas. Also F is the radiative flux, and k is the 
opacity, B is magnetic field and J is the current density, 
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i.e. J = V x B. We also assume that n is constant, and 
for simplicity, the diffusion of the radiation is neglected. 



p'g 



1 



(16) 



3 Linear Perturbations 

For doing a linear analysis, one should first specify the 
unperturbed properties of the system. We suppose that 
the streaming takes place in the a;— direction with a 
velocity U(z) as 



e r = 3p' r . 



(17) 



Having equations ©, ([16]) and ([I?]), then the lin- 
earized energy equation (fT2")l simplifies to 



p'r 



(18) 



U(z) 



U , 
-U , 



z>0 
z < 



where Uq is constant and all the unperturbed physical 
quantities are assumed spatially uniform in both the 
upper and the lower layers and we assume the initial 
magnetic field is along the streaming, i.e. B = Bgi 
where Bo is constant. In our configuration, the transi- 
tion layer at z = has a negligible thickness. Now, we 
can linearize the basic MHD equations. 

We perturb the physical variables as C (x, z, t) = Co + 
C,'{x,z,t), where Co is the unperturbed quantity and 
ICI "C | Co |- We substitute the perturbed variables as 



C(x,z,t) = az)e Hujt+k * x \ 



(8) 



where to is the growth rate and k x is the wavenumbcr 
of the perturbations. Thus, for a magnetized stream- 
ing fluid the basic equations (JTJ), (0), ([3]), g]), ^ arc 
linearized as 



du' 

iipp' + p (ik^u' x + —1) 
dz 



iippou^ = -ik*(p l+Pr), 



0, 



(9) 



(10) 



, , dp' dp' r , B , dB' . . , „ , 



dz dz Att dz 



du' 



itpE 1 + (E + P )(ik x u' x + ~r) = 0, 

dz 



(12) 



where 

£=ki+40), 



(19) 



and /3 — por/pog and we assume that the radiation 
pressure is smaller than the gas pressure which implies 
0<j8<l and !<£<§. 

Considering equations (fT0|) and (fT5|) . we can also 
rewrite equation ([9]) as 



-Pot 



du' 



h.2 

h 



^ (1+0 



(20) 



Substituting the above equation and equation (|18[) and 
equations (|X3j) and (fT4| into equation (fTT|) , the following 
differential equation for u' z is obtained 



dz 2 
where 



0. 



(21) 



1- 



4irtp 2 p 



cg(i + 0' 



47r<^ 2 po' 



(22) 



This is the main equation for analyzing KH instabil- 
ity, in which both the magnetic field and the radiation 
arc considered. The general solutions for the upper and 
the lower layers are written as 



u' lz = C ie - qiz + C[e +qu 



*2z 



C 2 e +92Z + C' 2 t 



z > 
z < 



UpB' x = -B 



dz ' 



iipB' z = -iB k x u' z , 

where ip = u> + k x Ua. Also, we have 



c 2 J, 



(13) 



(14) 



(15) 



where C\, C 2 , C[ and C' 2 are constants to be deter- 
mined based on the boundary conditions. Here, the 
parameters q\ and q 2 are corresponding to the upper 
and the lower regions, respectively. Without losing the 
generality of the problem, we assume both these param- 
eters have positive real parts. Once we obtain unstable 
modes, it will be checked if our assumption regarding 
to qi and q 2 are satisfied. 
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One boundary condition is that the flow should not 
diverge at infinity, i.e. u' z = when z —> ±00. Thus, 



c ie - qiz 

C 2 e +q2Z 



z>0 
z < 



Another boundary condition is the continuity of the 
displacement across the interface z = 0. So, we obtain 



Ci 



C2 

¥>2 



(23) 



Finally, the third boundary condition is the continu- 
ity of the total pressure at the interface. Now, we can 
obtain the dispersion relation. We assume the initial 
density and the sound speed at the upper layer are p\ 
and Cat, respectively. Also, p 2 and c s2 are the initial 
density and the sound speed at the lower layer. There- 
fore, the dispersion equation is obtained 



a(.x + l) 2 [l- 



1 



(x-l) 2 [l 



M| B (x-l)2 J 



1 - Mfezll! Ml B ( x -iy 



(24) 



where M\ and M 2 are Mach numbers in each layer, i.e. 
All = Uo/csi and M 2 = Uo/c s2 . Also, x is the dimen- 
sionless growth rate which is defined as x = co/h^Uo 
and the density contrast of the layers is denoted by 
a = p\j pi- Also, magnetic Mach number are defined 
as M\b = Uq/via and M 2 b = Uq/v 2 a, where via and 
v 2 a are the Afven speed in each layer. If we expand 
the above equation, then a six order polynomial is ob- 
tained. We solve this equation numerically to find the 
unstable modes. 

As another illustrative case, we consider two stream- 
ing layers where in one medium the radiation field is 
negligible but in the another medium the dynamical 
role of the radiation is considered. If we don't have ra- 
diation in fluid 2, the final dispersion relation becomes 



1 



1 



a(x + l) 2 [l 
1 



1 A/?0H-i)' M?Jx + l) 2 
(x-l) 2 [l 



MUx-l)^ 



i-M$(x-iy M$ B (x-iy 



= 



(25) 



We can also simplify the above relations by neglect- 
ing the magnetic field. When radiation field exists in 
both non-magnetized layers, one can simplify equation 
(121 as 



a(x + l) 2 



I) 2 



1+S 



0. 



A/j(z-l) 2 



(26) 



Also, for the non-magnetized case equation (|25[) reduces 
to the following equation, 



a{x + l) 2 



l _ A/?(x+i) 2 ^1 _ M'i{x - If 







(27) 



Equations ([24]), (|25j), (J26J) and (J27J) are our dispersion 
relations for analyzing KH instability under different 
conditions. In the next section, we study the unstable 
perturbations numerically. 



4 Analysis 
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Fig. 1 Growth rate of the unstable perturbation versus 
Mach number for different values of the parameter £ when 
Mi = 0.7 and a = 2. Each curve is labeled by the corre- 
sponding value of ^ and the dashed-dot curve shows growth 
rate for a case without radiation. 

First, we consider a non- magnetized system but with 
a radiation field. In this case, only one root of equa- 
tion (f2l)]) is acceptable and corresponds to the instabil- 
ity. Figure [1] shows growth rate of the unstable per- 
turbation versus Mach number M 2 for different values 
of the radiation parameter £. Here, we have Mi = 0.7 
and a = 2 and the dashed-dot curve corresponds to a 
case without radiation. As we see the curve without 
the radiation field has the lowest growth rate and with 
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Fig. 2 Growth rate of the unstable perturbation versus 
Mach number for different values of the parameter f for 
different magnetic Mach numbers, i.e. M\b = 2 (upper), 
Mis = 3 (middle) and Mis = 4 (bottom). Other input 
parameters are Mi = 0.7 and a = 2. Each curve is labeled 
by the corresponding value of £ and the dashed-dot curve 
shows growth rate for a case without magnetic field and 
with the radiation. Also, dot curves show growth rate for a 
magnetized case but without the radiation. 
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Fig. 3 Growth rate of the unstable perturbation versus 
Mach number for different values of the density contrast a. 
Each curve is labeled by the corresponding value of a and 
for three values for the radiation parameter, i.e. £ = 0.4 
(upper), £ = f.O (middle) and £ = 1.6 (bottom). Dashed- 
dot curve shows growth rate for a case without the radiation. 
All plots are corresponding to Mi =0.7 and a — 2. 



increasing the strength of radiation field, the system 
becomes more unstable. In other words, dynamically 
important radiation field has a destabilizing effect on 
the growth of the unstable KH perturbations. How- 
ever, the effect of the radiation is more significant at 
the larger Mach numbers. 

In Figure [2J we explore the combined effects of the 
magnetic field and the radiation. In order to make eas- 
ier comparison, growth rate corresponding to a mag- 
netized case without radiation is shown by dot curve. 
Top, middle and bottom plots arc corresponding to 
M\b = 2,3,4, respectively. While radiation pressure 
tends to destabilize the system, we see that magnetic 
field has a stabilizing role. All growth rate curves shift 
downward because of the magnetic field. In particular, 
when the magnetic Mach number M\b is larger mag- 
netic field tends to compensate the destabilizing role of 
the radiation field. 

In Figure [3l the effect of the density contrast a on 
the unstable KH perturbations is shown. In each plot, 
growth rates corresponding to the both cases with and 
without radiation are shown. The radiation parameter 
£ has three values, i.e. £ = 0.4 (top), 1.0 (middle) and 
1.6 (bottom). Each curve is labeled by the correspond- 
ing value of a. For a given radiation field, as the density 
contrast decreases, the effect of the radiation becomes 
more significant. In particular, for our input parame- 
ters, when the density contrast is a = 6 the difference 
between the cases with and without radiation is negligi- 
ble unless the Mach number becomes large where again 
the destabilizing role of the radiation field appears. 

5 Discussions 

We studied compressible KH instability in the presence 
of a dynamically important radiation field and the mag- 
netic field. Irrespective of the existence of the magnetic 
field, our results show that radiation field has a desta- 
bilizing effect on the unstable KH perturbations. One 
can simply define an effective pressure by adding ther- 
mal and radiation pressures. This will give us an effec- 
tive sound speed. Thus, in the presence of the radiation 
field the effective sound speed increases comparing to a 
case without radiation which then implies a less com- 
pressible system. In other words, dynamical role of the 
radiation field leads to a reduction of the compressibil- 
ity of the system, and so, KH perturbations have larger 
growth rates. Obviously, we do not expect to see the dy- 
namical effect of the radiation field on the incompress- 
ible KH perturbations. On the other hand, magnetic 
field has a stabilizing role. Thus, as we showed, there 
is always a competition between the radiation field and 



the magnetic field to control growth rate of the unsta- 
ble KH perturbations. We also found when the density 
contrast of the layers is smaller, the dynamical role of 
the radiation is stronger. But as the density contrast 
increases, the effect of the radiation field becomes more 
significant at the larger Mach numbers. 

Our analysis is restricted to the linear regime and 
so one can not apply the present results to the real 
systems to obtain the mass and energy exchanges at the 
interface. But we expect stronger turbulence via KH 
instability at the nonlinear regime when dynamical role 
of the radiation becomes significant. One of the sources 
of the observed turbulence in ISM is the injection of 
mass and energy from the stellar outflows and jets. KH 
mechanism provides a way for these exchanges which 
would be amplified in the presence of the radiation field. 
Also, it is indeed likely that radiation pressure could 
have a dynamical effect in the boundaries of radiation 
driven bubbles where KH instability may occur. As 
a next step for the future work, numerical simulation 
of KH instability in the presence of the radiation is 
required to address these issues. 
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